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Abstract 

In this paper, we report the results from the analysis of a numerical 
model used for the design of a magnetic linear actuator with applications 
to fin-based locomotion. Most of the current robotic fish generate bending 
motion using rotary motors which implies at least one mechanical conver- 
sion of the motion. We seek a solution that directly bends the fin and, 
at the same time, is able to exploit the magneto-mechanical properties 
of the fin material. This strong fin-actuator coupling blends the actuator 
and the body of the robot, allowing cross optimization of the system's 
elements. 

We study a simplified model of an elastic element, a spring- mass sys- 
tem representing a flexible fin, subjected to nonlinear forcing, emulating 
magnetic interaction. The dynamics of the system is studied under un- 
forced and periodic forcing conditions. The analysis is focused on the limit 
cycles present in the system, which allows the periodic bending of the fin 
and the generation of thrust. The frequency, maximum amplitude and 
center of the periodic orbits (offset of the bending) depend directly on the 
stiffness of the fin and the intensity of the forcing; we use this dependency 
to sketch a simple parameter controller. Although the model is strongly 
simplified, it provides means to estimate first values of the parameters for 
this kind of actuator and it is useful to evaluate the feasibility of minimal 
actuation control of such systems. 

1 Introduction 

In the last two decades underwater fin-based propulsion has been a topic 
of intense research. Theoretical models and simulations based on a mul- 
tiplicity of numerical methods have been developed and reported 
9 ( If, 12, 15, 18 as well as experiments performed on artificial platforms 
and with animals [T1Ffl ll0l[T3] . Theoretical studies such as [6l[l5] focus 
on models of the viscoelastic body, while [Hl l 1 1IIT5] treat the body-fluid 
coupling and the emergence of locomotion. In 2 ,9, 18 the main interest is 
to understand how passive thrust is generated in vortex wakes, a question 
that remains open. Despite the considerable work done, experimental re- 
sults are not fully understood and several numerical simulations are yet to 
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be validated. Influenced by the activity on the field, we address a parallel 
problem related to robotics, actuation. 

Fins as a tool for locomotion offer several appealing properties with 
respect to propellers. From an environmental point of view fins reduce 
sound pollution characteristic of propellers [16] • Additionally, a flexible 
body offers the possibility of extracting energy from the environment, as 
shown in the technological study presented in [4]. This idea is related 
to the fact that trouts exploit vortex wakes (shed by obstacles in the 
flow) to reduce the cost of swimming [5] (see [T3] for a review). These 
aspects are of primary relevance in situations where low environmental 
impact and mimicry are important, as in pipes maintenance routine, or 
for underwater life observation (scientific naiveness may make us forget 
about military applications, we provide this short caution.). Moreover, 
in environments where moving parts may be clogged up due to fouling, 
rotatory propellers may be unfeasible for locomotion. 

In the field of bio-inspired robotics, actuation for swimming robots 
using fin-like propulsion is usually implemented by the use of rotatory 
electric motors to control the fin swing angle. However, The presence 
of motors hinders miniaturization and integration of actuators into the 
robot structure, a requirement for flexible machines with deformable body. 
At the same time this imposes a mechanical conversion of rotations into 
oscillatory linear motion, the complexity of which becomes an issue at 
smaller size. Moreover, minimizing or eliminating the number of moving 
parts required for actuation simplifies maintenance to a large extent. To 
tackle the mentioned difficulties we are directing our research towards new 
ways of actuation. Herein we report a simple mathematical model and 
the numerical analysis of a possible alternative. 

We start by considering a robot composed of a hull and a fin attached 
to it. The fin is modeled as an elastic beam (see [TS] or [5] for more detailed 
models), which we want to set into oscillatory motion. In the setup shown 
in Fig. [T] we choose to support the beam at two points. The first support 
is at one edge of the beam and stands for the hull of the robot. The second 
support is placed at some intermediate point in the beam. The section 
of the beam beyond this second support is meant to generate the thrust 
by interacting with the surrounding fluid. The actuation is done in the 
section of the beam between the two supports by means of a combination 
of permanent magnets (one of them attached to the beam) and solenoids. 
In the configuration chosen, the permanent magnets serve to increase the 
compliance of the system and to reduce the force that needs to be actively 
applied by the solenoid. The distance between the supports defines the 
rigidity of the actuated section and could be tuned for optimal energy 
transfer. Similar working principles are described in patents of electric 
razors, and of active dampers of oscillations for digital cameras lenses 
(in these contexts the actuator is often called motor or electromagnetic 
spring). Similarly, the control of the resonant modes of a structure is 
a commonplace problem in structural dynamics [19] ■ It is noteworthy 
that all these techniques exploit (to be more efficient) or require (to be 
implemented) knowledge of the resonant modes of the system under study. 
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Figure 1: Schematic of the system described by equations (P). A beam is used 
to model the fin and it is simply supported at two points. The actuation is done 
by means of a combination of permanent magnets and solenoids. The distance 
between the supports defines the rigidity of the actuated section. 



2 Dynamic Model 

The displacement of the magnet in the fin (we will refer to this point as 
the fin magnet), can be modeled by a spring-mass system under the effect 
of an external force field. Considering only one dimensional motion, the 
system is written as, 

x = v 

. F 2 (x) + Fi(a;) + F sl {x) + F s2 {x) K T C 1 ) 

V = — — — — X V, 

m mm 

where x is the displacement of the fin magnet, m is an effective moving 
mass, K represents an effective elastic constant of the fin setup and Y is 
used to include dissipation. The Fi and F 3 i terms are forces acting on the 
fin due to the external magnets and the solenoids, respectively. 

Magnetic forces can be highly complex; to keep our model as simple as 
possible we approximate each magnet as a point magnetic dipole, which 
is a good approximation when the distances are significantly bigger than 
the size of the magnet in the direction of the magnetization [20] . In this 
situation the force can be expressed as follows, 

Ci 

Fi{x) = —. l —-5;sign(x-Xi), (2) 

(a; — Xi) 

where d is a constant that depends on the magnetic moments of the 
magnets and their geometry, positive values represent attractive forces and 
negative values repulsive forces. In the case of the solenoids, this constant 
depend also on the current, i.e. C 3 i(I) = C 3 il(t) (the index s refers to 
solenoid). The position of the external magnet (or solenoid) measured 
from the rest position of fin magnet is Xi. Henceforth we define k — K /m, 
7 = r /m, Ci = c i/m and c 3 i — c a»/m. Additionally, we assume that the 
deflected fin does not reach the external magnets, in mathematical terms 
this is expressed as i 6 (xi,X2)- 

Next we study the dynamics of the system without actuation, I(t) = 0. 
The expression for the fixed points x* is obtained by equating system fTJ 
to zero. The second equation yields 

C2 (a;* — x{) a — ci (x* — 22)" — kx* [(x* — xi) (x* — 12)]° = 0, (3) 
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where the assumption x € (xi, X2) was used to determine the signs. 



Linear Stability Analysis. To classify the fixed points, we calculate 
the trace and determinant of the 2x2 Jacobian matrix J of ([1]). These 
are given by 

Tr(J) = -7 (4) 

C2 Ci 



Det(J) = k + a 



[X — X2) {X — XX 1 



(5) 



In general the fixed points of the system will be saddle-nodes, centers 
or spirals, depending on the value of the parameters 7, k, a and Xi. 
However, the position of the fixed point (i.e. the solutions of Q) are 
independent of 7. 

To proceed with the analysis we introduce further assumptions. The 
exponent a depends on the arrangement of magnets [20] . Here we will 
consider identical cylindrical magnets placed symmetrically with respect 
to the rest position of the fin magnet and with dipoles parallel to it (at- 
tracting); hence a — 4, c\ = C2 = c > and X2 = —xi — xq > 0. By 
neglecting dissipation, i.e. 7 = 0, we set the trace of the Jacobian to zero. 
Consequently, the fixed points are either saddle-nodes or centers, depend- 
ing on the sign of ([5}- Using these assumptions to simplify the equality 
@ we obtain, 



8cxn 



+ xl) - k (x* 2 - xq) 4 ] = 0, (6) 



rendering evident that x* — is one of the fixed points, in consequence 
of the symmetry of the problem. The determinant ([5]) at this point is, 

Det(J)|«.= = A-84. (7) 
x 

which is positive for c /k < x o/&, and the origin is a center. Although any 
real system will not show centers without actuation (due to dissipation) 
their position will match the pole of the spirals observed. 

It can be shown that the nonzero solutions of (|6j) are saddle-nodes. 
Displacements beyond the saddle-nodes will bring the fin magnet into a 
region where the attraction is stronger than the elastic restitution, causing 
the fin to stick to the closest magnet. The saddles establish a natural limit 
for the maximal amplitude of the orbits of the system. To illustrate these 
ideas, we show in Fig. [2] three plots of the polynomial defined by ((SJl 
for different values of the ratio c /k, together with phase portraits of the 
system. The figure depicts the trade-off between the rigidity of the fin and 
the interaction of the fin magnet and the permanent magnets. Keeping 
the xo fixed, the stronger the magnets (or the more compliant the fin), 
the smaller the region where the system can present stable orbits. At 
the critical ratio c /k = x o/s, the saddle-nodes collide at the origin and the 
center is transformed into a saddle-node. 



3 Actuation and Control 

As mentioned before, dissipation will reduce the amplitude of the oscilla- 
tions. Therefore, to keep the system close to the desired trajectory we need 
to pump energy into it. To do this, we have placed a solenoid surrounding 
both magnets such that they can increase or decrease the interaction with 
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Figure 2: Plot of the polynomial denned in (J6]) for different values of c /k. The 
star symbols mark the position of the fixed points. The phase portraits to 
the right show that the saddle-nodes define a limit for the amplitude of the 
orbits. The figure illustrates the trade-off between the rigidity of the fin and 
the intensity of the magnetic interaction. 



the fin magnet. Both solenoids are constructed similarly but arranged 
anti-parallel to each other (c S 2 = — c s i = c s > 0). We place them as close 
as possible to the fin magnet, for example near the saddle-nodes of the 
system. In order to determine the parameters of the solenoid needed to 
drive and control the system, we use a simple PID controller which can 
regulate the applied force. To this end, we rearrange the terms of (JTJ) and 
write them as, 



- + F c (t) — kx — jv 



(x - x ) (x + x ) 
F c {t) = k p e(t) + k d e(t) + h I e(s)d. 



Is 



(8) 



o 



e(t) = x d (t)-x(t). 



Where Xd(t) is the desired displacement. The effort required to drive the 
system into steady oscillations, depends on the appropriate choice of k p , 

kd> ki. 



4 Parameter values 

In the following sections we define the values of the parameters used for 
the numerical simulation of the platform. We give a brief description on 
the assumptions and criteria used to select them. In Table[T]we summarize 
the information. 

Magnets. The constant for the force (J2J) in the case of cylindrical mag- 
nets magnetized along their length £ is, 

where fj,o is the permeability of vacuum, R is the radius of the magnet and 
B r is the remanence of the magnet (a value available from manufacturers). 
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The factors in parenthesis represent the magnetization and the volume of 
the magnet, they are squared because we are assuming both interacting 
magnets are equal. Using values from commercially available Neodymium 
magnets we have calculated C = 2.460 x 10~ 10 N ■ m 4 . 



Elastic constant. As can be seen in Fig. [T]we model the fin using 
a beam simply supported in two points. The supports are separated by 
a distance L. The fin magnet is placed at the point y and the elastic 
stiffness there can be written as, 

Q T 

where E is the Young's modulus of the material and J is the second 
moment of area of the beam. As discussed before, the behavior of the 
system depends on the relation between the elastic constant of the fin 
and the strength of the magnets. The current setup allows tuning the 
elasticity of the fin by setting different materials and profiles of the fin, 
or by moving the fin magnet, or by changing the distance between the 
supports. If needed, the setup could be transformed into a cantilever by 
removing the second support. To provide good ranges of elasticity we use 
Low-density polyethylene plastic with E ~ 0.2 x 10 9 Pa and L = 30mm. 
The fin has width and thickness of 10mm x 0.5mm, respectively. The 
fin magnet is placed in the middle of the two supports. This values yield 
K = 37 \03N ■ m" 1 . 



Damping and mass. When a body moves in a liquid, it transfers ki- 
netic energy to the surrounding fluid reducing the acceleration it presents 
corresponds to the one observed on a body with higher mass. This phe- 
nomena, known as added mass, can be estimated using models as the one 
presented in |21| . However, we postpone a detailed description for future 
work and simply consider a total mass 300 times bigger than the mass of 
the fin and the magnet together, m = 84. 3g. Estimation of the damping 
r without an experimental setup is not straightforward, therefore we use 
damping ratios in the range Q € [0, 0.5]. Where Q = means no damping 
and a value Q = 1 corresponds to critical damping. 

The table below summarizes the value of the parameters. Though the 
values are reasonable, we do not expect them to correspond to any real 
device and corrected ones will come from a future validation process. 



Parameter 


Value 


Units 


C ( C /m) 
k (IC/m) 
d (V m ) 
x 


2.919 x 10" y 
439.3 
[0, 20.96] 

i x nr 2 


N ■ to 4 • kg^ 1 
N ■ mT 1 ■ kg^ 1 
N ■ s ■ tn~ ■ kg^ 1 
m 



Table 1: The table summarizes the values of the parameters used for the nu- 
merical results reported in the text. Values for d correspond to the range of 
damping ratio Q 6 [0,0.5]. 
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5 Results 



Phase space and time series. Numerical results for the undamped 
system are presented in Fig. We take three initial conditions on the 
region of the phase space to study. All the initial conditions start with 
zero velocity, i.e. they lay on the horizontal axis. It is important to 
note that the time series of the fin displacement clearly show different 
frequencies. This is due to the attraction of the magnets, the higher the 
initial displacement the lower the frequency of the orbit. These results 
are shown in detail in Fig. [3] For each initial condition we plot the power 
spectrum of the signal and it is visible how the main component decreases 
at higher amplitudes. 



Position [mm] 




Figure 3: Trajectories in phase space and time series for the undamped sys- 
tem starting from three different initial condition. The frequency of the signal 
decreases with the amplitude due to the interaction of the magnets. 

The offset of the oscillations corresponds to the position of the center. 
By breaking the symmetry of the system, either by setting ci 7^ C2 or 
by feeding constant current to the solenoids, we can move the center off 
the origin. This could be required for turning maneuvers or useful for 
initiating oscillations. In Fig. [5] we show how the center and the saddles 
move for different values of ci and ci (or increasing current). 

Actuated system. The controller is provided with reference signal 
of the form Acos(wt + <f>). The amplitude A is taken inside the bounds 
defined by the saddle-nodes. Given a value for the amplitude, we use the 
curve shown in the inset of Fig. [4] as a frequency lookup table. The initial 
phase (j> is calculated from the initial conditions. We took the same initial 
conditions as in the previous section, namely x(0) = A, 0.5 ■ A, 1.25 • A. 
However simple, the control technique shows a remarkable performance as 
can be appreciated in Fig. [U In the figure the reference trajectory is the 
dashed line, for the case without damping, the actual trajectories of the 
system overlap with. On the left panels the time series of the displacement 
are shown and compared with the reference amplitude. 



Solenoid. If we take the maximum of the output from the controller, 
F c (t), and limit the maximum current fed to the solenoids to / < 20mA, 
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Figure 4: Variation of the natural frequency with the amplitude of the oscilla- 
tions. The power spectrum of the orbits is plotted, the behavior of the main 
component is shown in detail in the inset. 
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Figure 5: Control of the position of the center and saddles. The center moves 
symmetrically around the origin for differences Ac between the magnetic con- 
stants. The maximum amplitude is also compromised, because the saddle on 
the side of the stronger magnet come closer to the center. 



we can find a suitable expression for calculating the ideal number of turns 
N of a coil. The net force generated on the fin magnet is given by the 
sum of forces due to the 2 coils, the wiring of which ensures that the 
force generated by each on the fin magnet is additive. This force can be 
expressed as, 



F m = g {B r R 2 e) (NIA) G(x r 



G(x m ) = 



(x m + x ) 4 + (x m - x ) 4 



(11) 



2\4 
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Figure 6: Trajectory in phase space of the actuated system for three different 
initial conditions close to the desired amplitude. On Top the results without 
damping and below them, with damping ratio of Q = 0.5, which gives an 
unfeasible solenoid. In the second case the reference trajectory in phase space 
is shown in dashed line. On the left we show the time series of the displacement 
of the fin. 

where B r ,R and £ are as in [5] and x m is the displacement corresponding 
to the maximum force F m . Solving for N and replacing with the corre- 
sponding values we get N £ [600, 1000], for Q G [0,0.2]. Higher values of 
damping impose too many turns on the solenoid. Although the model is 
not yet validated this result is encouraging. 

6 Discussion and conclusion 

In this short report we have presented the first step towards the design and 
construction of a novel actuator for small swimming machines. Though 
we have used rough models, the results show that we are pointing in 
the right direction. The design shown here is not necessarily the best in 
reducing the actuation needed. For example, one could think of using 
the instability of the center to initiate motion, by forcing the system 
through its bifurcation. This could be achieved by on-line modification of 
the distance between the supports or by bringing the permanent magnets 
closer. Additionally, placing a permanent magnet perpendicular to the 
plane at the origin could be used to further reduce the frequency of the 
orbits or to control the offset in a more sensitive way than the one shown 
here. 

Our model includes dissipation proportional to the velocity and there- 
fore the role of dissipation is marginal. More detailed models of the fluid 
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dynamics and the bending of the fin will surely bring dissipation into a 
more primary role in the behavior of the system. In addition, thrust, heat 
dissipation and energy consumption could be estimated in such multi- 
physics models. 

We have shown how a simple PID controller could perform reasonably 
when information about the phase portrait of the system, like the depen- 
dence of frequency with amplitude, is included. The use of adaptable fre- 
quency oscillators [T7] or standard model-based controllers (like feedback 
linearization), could improve performance and reduce these requirements. 
Additionally, a system that is too flexible does not possess orbits without 
a controller. Such a controller would requires large amount of actuation, 
since it is forcing the system to behave unnaturally. Therefore existence 
of orbits can be exploited to reduce energy consumption. This stresses the 
fact that passive dynamics are a key to improve the way we control and 
design our robots. Controlling the force between solenoids and moving 
magnets, brings several challenges on the design of the electrical circuits 
due to the changes in impedance. Another interesting aspect of the prob- 
lem that will be addressed in further studies. We understand that results 
obtained solely from simulations are as 'words without actions', however 
the use of simple models can help us evaluate the feasibility of certain 
designs. In our particular case, a device with low friction could be driven 
with a tuned PID, a frequency lookup table and a solenoid with 800 
turns, consuming about 20mA. 
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